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Abstract 

We develop and assess the potential of several powerful techniques, designed to 
investigate the details of reionization. First, we present a procedure to probe the 
neutral fraction, xm, using the Lyman alpha transmission statistics of high-redshift 
(z ^ 6) sources. We find that only tens of bright quasar spectra could distinguish 
between xhi ~ 1 and xhi ^ 10~ 2 . A rudimentary application of such a technique on 
quasar SDSS J1030+0524 has yielded compelling evidence of a large neutral fraction 
(xm ^ 0.2) at z ~ 6. We also generate the observable, high-z supernovae (SNe) rates 
and quantify the prospects of detecting the suppression of star-formation in low- 
mass galaxies at reionization from such SNe rates, specifically from those obtainable 
from the James Webb Space Telescope (J WSJ). Our analysis suggests that searches 
for SNe could yield thousands of SNe per unit redshift at z ~ 6, and be a valuable 
tool at studying reionization features and feedback effects out to z ^ 13. 
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1 Introduction 



The epoch of cosmological reionization is a significant milestone in the history 
of structure formation. Despite recent observational break-throughs, the de- 
tails of the reionization history remain poorly determined. The Sloan Digital 
Sky Survey (SDSS) has detected large regions with no observable flux in the 
spectra of several z ~ 6 quasars (e.g (|l|)). The presence of these Gunn- Peterson 
(GP) troughs set a lower limit on the volume weighted hydrogen neutral frac- 
tion of xhi > 10~ 3 (0), implying a rapid evolution in the ionizing background 
from z = 5.5 to z ~ 6 (e.g. (|l|)), and suggesting that we are witnessing the 
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end of the reionization epoch, with the IGM becoming close to fully neutral 
at z ~ 7. On the other hand, recent results from the Wilkinson Microwave 
Anisotropy Probe (WMAP) have uncovered evidence for a large optical depth 
to electron scattering, r e ~ 0.17 ± 0.04 (2) in the cosmic microwave back- 
ground anisotropics. This result suggests that the universe was already highly 
ionized at redshifts as high as z ~ 15 - 20. Various feedback mechanisms have 
been proposed to regulate the evolution of the ionization state of the IGM, 
with no clear consensus on a favored plausible physical model (e.g. (3)). It 
is evident that more tools designed to probe the reionization epoch would 
be of great value in pinning-down the reionization history and increasing our 
understanding of the early universe. 

The rest of this contribution is organized as follows. In § 2 we discuss how the 
transmission statistics of high-redshift spectra can be used to place constraints 
on the neutral fraction and size of HII regions. In § 3, we describe how the 
extended dynamical range provided by the Lyman a and Lyman f3 absorption 
allowed us to model the gross observed features of quasar SDSS J1030+0524. 
In § 4, we predict observable high-redshift SNe rates and asses their usefulness 
in the detection of reionization features. Finally, in § 5, we offer our conclu- 
sions. Most numerical estimates presented here assume standard concordance 
cosmological parameters, (fl\, f^M, n , cr 8 , H ) = (0.73, 0.27, 0.044, 1, 0.9, 
71 km s _1 Mpc -1 ), consistent with the recent WMAP measurements. Unless 
stated otherwise, all lengths are quoted in comoving units. 



2 Mock Spectral Analysis 

In principle, the quasar's absorption spectrum contains a full record of the 
neutral fraction as a function of position along the line of sight. Since high- 
redshift sources sit in their own highly ionized Stromgren spheres, the total 
Lyman a optical depth at a given observed wavelength, A Q bs, can be written 
as the sum of contributions from inside and outside the Stromgren sphere, 
T hya = t r + t d- These two contributions are shown in Figure 2, originating 
from a hydrodynamical simulation (further details of the simulation can be 
found in (5)). The residual neutral hydrogen inside the Stromgren sphere 
resonantly attenuates the quasar's flux at wavelengths around X a (l+z), where 
X a = 1215.67 A is the rest-frame wavelength of the Lyman a line center. As a 
result, t r is a fluctuating function of wavelength (solid curve), reflecting the 
density fluctuations in the surrounding gas. In contrast, the damping wing of 
the absorption, td, is a relatively smooth function (dashed curve), because its 
value is averaged over many density fluctuations. The damping wing optical 
depth is a strong function of the size of the Stromgren sphere, Rs, and the 
hydrogen neutral fraction outside the Stromgren sphere, xhi- 
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Although there is a wealth of information in these two components, it is sys- 
tematically challenging to separate them and hence extract relevant parame- 
ters. In we study the feasibility of statistically extracting such parameters 
from high redshift spectra. The free parameters in our analysis are Rs and 
xhi- Note that changing Rs moves the dashed (td) curve left and right, while 
changing % moves this curve up and down in Figure 2. 

The analysis can be summarized as the following. We start with a simulated 
observed spectrum, generated using a randomly chosen line of sight (LOS) 
from a cosmological hydrodynamic simulation. Then we guess values for the 
radius of the Stromgren sphere, R' s , and the IGM hydrogen neutral fraction, 
x' m . Next we approximate the amplitude of the source's intrinsic emission, A', 
implied by the choices of R' s and x' m , using the red side of the Lyman a line 
where resonance absorption can be neglected. From the observed spectrum, 
we divide out the assumed intrinsic emission (A f x known spectral shape), and 
the assumed damping wing flux decrement, e TD ^ obs ' R s> x ra> , calling the result 
S'(\ hs)- If our choices of R' s and x' m were correct, S"(A bs) should represent the 
resonance absorption flux decrement alone. Hence, we compare a histogram of 
the implied resonance optical depths, — ln[S"(A b s )], to the known histogram 
of resonance optical depth (obtained from the simulation). We then repeat 
this procedure with different choices of R' s and x' m , finding the ones whose 
implied resonance optical depths most closely match the known histogram. 

A few of the resulting histograms are shown in Figure 1. The template distri- 
bution of resonance optical depths obtained from our simulation is shown in 
the top left panel. Distributions derived from our inversion analysis explained 
above are shown in the other panels, with the correct parameter choice in 
the top right and incorrect choices in the bottom two panels. We test the hy- 
pothesis that the top right, bottom left and bottom right histograms (among 
many others in parameter space) were drawn from the distribution in the top 
left. We find that the top right panel is consistent with being drawn from the 
distribution in the top left, and that the bottom panels are not. 

Additionally, we have explicitly incorporated into our analysis an error in the 
intrinsic emission template, consisting of either an uncertainty in its spectral 
power-law index, or Gaussian, uncorrelated, pixel-to-pixel variations at each 
wavelength. With both of these errors, we find that a neutral universe can 
be statistically distinguished from a xhi = 0.008 universe in our parameter 
space, using tens of bright quasars, a sample that can be expected by the 
completion of the Sloan Digital Sky Survey. Alternatively, similar statistical 
constraints can be derived from the spectra of several hundred sources that 
are ~ 100 times fainter. For example, the Large-aperture Synoptic Survey 
Telescope (LSST) should be able to deliver many new faint quasars that could 
serve as targets for low-resolution spectroscopy. Furthermore, if the size of 
the source's Stromgren sphere can be independently constrained to within ~ 
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Fig. 1. Histogram of the template distribution tr (top left), and histograms of 
the derived t' r distribution: (xm,Rs) = (1,43 Mpc) (correct values) (top right), 
(xm,R s ) = (0.2,43 Mpc) (bottom left), (x m ,Rs) = (1,85 Mpc) (bottom right). We 
test the hypothesis that the top right, bottom left and bottom right histograms 
(among many others in parameter space) were drawn from the distribution in the 
top left. We find that the top right panel is consistent with being drawn from the 
distribution in the top left, and that the bottom panels are not. 

10% (such as with the method presented in (0)), the analysis presented here 
can distinguish between sources embedded in an IGM with 10~ 3 < xui < 1, 
using a single source. We plan to perform such analysis on the current sample 
of high-redshift sources. 



3 Case of SDSS J1030+0524 



The challenges discussed above in extracting relevant parameters from high- 
redshift spectra are partially aggravated by the limited dynamical range probed 
by observations. For example, the sharp rise in r D at wavelengths A Q b s ^ 8720 
A in Fig. 2 is a unique feature of the boundary of the HII region, and cor- 
responds to absorption of photons redshifting into resonance outside of the 
Stromgren sphere. The detection of this feature has generally been regarded 
as challenging: since the quasar's flux is attenuated by a factor of exp(— r Lya ), 
an exceedingly large dynamical range is required in the corresponding flux 
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measurements. However, in (|6j) we show that simultaneously considering the 
measured absorption in two or more hydrogen Lyman lines can provide the 
dynamical range required to detect this feature. 

The observational input to our analysis is the deepest-exposure absorption 
spectrum of SDSS J1030+0524 available to date (@). The flux detection thresh- 
old in the Lyman a and Lyman (3 regions of this spectrum correspond to total 
Lyman a optical depths of Twi^a) ~ 6 and T lim ( Ly /3) ~ 11 respectively. To 
summarize the implied constraints, we have divided the spectrum into three 
different regions, shown in Figure 2. In Region 1, with A D b s > 8752.5 A, the 
detection of flux from SDSS J1030+0524 corresponds to the upper limit on 
the total optical depth TL ya < 6. Region 2, extending from 8725.8 A < A D bs < 
8752.5 A, is inside the Lyman a trough, but outside the Lyman f3 trough. 
Throughout this region, the data requires 6 ^ r Lya ^ 11. Region 3, with 
A bs < 8725.8 A has a lower limit r Lyce > 11. As defined, each of these 
three regions contains approximately eight pixels in the spectrum of SDSS 
J1030+0524. 

We modeled the absorption spectrum of SDSS J1030+0524, attempting to 
match these gross observed features. We extracted density and velocity infor- 
mation from 100 randomly chosen lines of sight (LOSs) through the simulation 
box (described in (0)). The density and velocity distributions are biased near 
the density peaks where a quasar may reside. However, the spectral region 
of interest lies well outside these biased regions, which, in the case of high- 
redshift quasars, extends only to ^ 1 proper Mpc, on average (jg). The use 
of randomly chosen LOSs is therefore an accurate statistical representation of 
the expected density and velocity fields. Along each LOS, we computed the 
Lyman a absorption as a function of observed wavelength. 

In this analysis, we had three free parameters: (i) the size of the ionized re- 
gion (Rs)'i (h) the fraction of neutral hydrogen outside it (xhi); (hi) the QSO's 
ionizing luminosity {L v ). To understand the significance of these three param- 
eters, note that changing R$ moves the dashed (td) curve left and right, while 
changing % moves this curve up and down in Figure 2; changing L u corre- 
sponds to moving the solid (r R ) curve up and down. Values of xhi < 1 imply 
the existence of a background flux of ionizing radiation, with an ionization 
rate of Ti2 ~ lO -5 ;^ 1 x 10 -12 s _1 . For consistency, we add this constant back- 
ground flux to that of the quasar (the latter dropping as r~ 2 with distance), 
somewhat reducing the neutral hydrogen fraction inside the ionized region. 

We evaluated tr and tq for each LOS in our simulation and for every combi- 
nation of Rs, Xm, and L u , attempting to match the gross features explained 
above. In order to be conservative, we allowed up to two pixels in each region 
to lie outside the allowed range of Lyman a optical depths; each LOS that 
had more than two pixels fail the above criteria in any of the three regions 
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was rejected. We find that a more stringent criterion of allowing only a single 
'faulty' pixel would strengthen our conclusions. 



The procedure outlined above turns out to provide tight constrains on all three 
of our free parameters simultaneously. The radius of the Stromgren sphere is 
limited to 42 Mpc < R,s < 47 Mpc, close to the previously estimated lower 
limits (jit lioh. We also find evidence that the IGM was significantly neutral 
at z ~ 6, with a ~ 1 a lower limit of xhi ^ 0.17. This result is derived from 
the observed sharpness of the boundary of the HII region alone, and relies 
only on the gross density fluctuation statistics from the numerical simulation. 
In particular, it does not rely on assumptions about the mechanism for the 
growth of the HII region. Finally, we find an emission rate of ionizing photons 
per second of (5.2 ± 2.5) x 10 56 photons sec -1 , which is between 2-10 times 
lower than expected (11; ll2T ). strengthening arguments that reionization at 



z ~ 6 is caused by the radiation from early stars, rather than bright quasars 
0). 

Our findings represent the first detection of the boundary of a cosmologi- 
cal HII region, and have several important implications. As was previously 
mentioned in § 2, tight constraints on the hydrogen neutral fraction can be 
extracted directly from the Lyman a absorption spectrum alone (p), but this 
requires tens of spectra to be statistically significant when R$ is not known. 
Incorporating an independent limit on Rs, such as the one obtained from this 
method, can reduce the number of required spectra to one. Tight constraints 
on xhi have also recently been obtained by adopting a proper size of Re ~ 4.5 
Mpc for SDSS J1030+0524, together with an estimate of its lifetime (Hi). Our 
direct determination of the Stromgren sphere size is only slightly larger than 
the previously adopted value, lending further credibility to this conclusion. 
The sharp boundary we detect also constrains the hardness of the ionizing 
spectrum of SDSS J1030+0524. We infer here a rise in the neutral fraction 
over a redshift interval of Az ~ 0.02, corresponding to a proper distance of 
~ 1.2 Mpc. The thickness of the Stromgren surface is approximately given by 
the mean free path of the typical ionizing photon. In order for this mean free 
path not to significantly exceed the observed thickness, the mean energy of 
the ionizing photon emitted by SDSS J1030+0524 must be < 230eV. This is 
a interesting limit, implying that the effective slope of the ionizing spectrum 
of this source above E = 13.6eV is softer than —d\n.L v /d\n.v = 1.07. 

In the future, given a larger sample of quasars at z > 6, it will be possible to 
use the method presented here to search for sharp features in the absorption 
spectrum from intervening HII regions, not associated with the background 
source itself. We plan study the utility of this approach in a future paper. 
For such external HII regions that happen to intersect the line of sight, both 
the red and the blue sides of each GP trough can, in principle, be detected, 
and can yield two separate measurements of the ionized fraction at different 
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Fig. 2. Model from a hydro dynamical simulation for the optical depth contributions 
from within (tr) and from outside (td) the local ionized region for a typical line 
of sight towards a z s = 6.28 quasar embedded in a fully neutral, smooth IGM, 
with Rs = 44 Mpc and /i on = 1- The dashed curve corresponds to td, and the 
solid curve corresponds to tr. The total Lyman a optical depth is the sum of these 
two contributions, TL yQ , = tr + td- The dashed-dotted lines demarcate the three 
wavelength regions used for our analysis described in the text. In our analysis, 
the optical depths are averaged over 3.5 A wavelength bins, which decreases the 
fluctuation of t#. For reference, the redshifted Lyman a wavelength is at 8852 A. 



points along the line of sight. The Universe must have gone through a tran- 
sition epoch when HII regions, driven into the IGM by quasars and galaxies, 
partially percolated and filled a significant fraction of the volume. The detec- 
tion of the associated sharp features in future quasar absorption spectra will 
provide a direct probe of the the 3D topology of ionized regions during this 
crucial transition epoch. Deep surveys equipped with sufficiently red niters 
(with instruments such as those being carried out by the VLT, and ultimately 
with deep surveys coverin g a significant portion of the sky, such as the survey 
proposed with the LSST; (|l5l)) will be able to deliver a large sample of bright 
z > 6 quasars needed for such studies. 
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4 High-Redshift SNe Rates 



During reionization, the ionizing background radiation heats the IGM and 
could dramatically suppress gas accretion and star formation in low-mass 
halos. As a result, reionization may be accompanied by a drop in the star- 
formation rate (SFR) and the corresponding supernovae rate (SNR). The size 
of the drop is uncertain, since the ability of the gas in low-mass halos to cool 
and self-shield against the ionizing radiation is poorly constrained at high 
redshifts. The redshift and width of the drop are also uncertain as the details 
of reionization are also unknown. In |l6h. we construct the observable, high-z 
SNe rates, and analyze the prospects of detecting a reionization feature with 
future SNe observations, specifically those obtainable with JWST. 

Such a drop in the SFR could be detected in the high-redshift extension of the 
'Madau' diagram (e.g. (fin) ), by directly counting faint galaxies (fill). However, 
these galaxies are dangerously close to the detection thresholds of forthcoming 
instruments, making a reionization feature detectable only in very optimistic 
scenarios. Alternatively, by analyzing the Lyman a absorption spectra of SDSS 
quasars, (IS) suggest that we may already have detected a drop in the SFR 
at z ~ 6, through the non-monotonic evolution of the mean IGM opacity. In 
order to improve on this current, low signal-to-noise result, many deep, high- 
resolution, high- z spectra of bright quasars would be required. Other events 
that may trace out the early SFR, such as gamma-ray bursts (GRBs) and SNe, 
could be better suited to detect the reionization drop, as they are bright, and 
could be used as tracers of the SFR in arbitrarily faint host galaxies, out to 
redshifts higher than the host galaxies themselves. Unfortunately, while GRBs 
are bright, they are rare. Based on the expected Swift GRB detection rates 
(e.g. Figure 6 in (jUJ), a drop in the GRB rate associated with reionization 
could be detected only at very low, ~ la, significance, and only if reionization 
occurred at z re ^ 7 - 10. On the other hand, SNe have the benefits of being 
both very bright (compared to galaxies at the very faint end of the luminosity 
function), and occurring much more frequently (compared to GRBs). 

There have been several previous studies of the expected early SNR. How- 
ever, these either did not focus on the expected rates at redshifts beyond 
reionization (e.g (|2lh). focused exclusively on very high-z rates from the first 
generation of metal-free stars (1221) . and/or did not address the impact of reion- 
ization on the SFR (e.g. (|23t \24\ ). Likewise, previous theoretical predictions 
of the high-redshift SFR (e.g. (|la)) assumed a fixed degree of suppression due 
to reionization, matching numerical simulations which did not include self- 
shielding (J25h. The main distinctions here are that we compute the expected 
SNe rates at a wide range of redshifts and we allow various degrees of pho- 
toionization heating feedback. We also supplement the standard halo-based 
estimates of the SFR with a more elaborate estimate of the corresponding 
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observable SNR, utilizing the observed properties of \ow-z SNe. 



We assume that the SFR at high redshift traces the formation rate of dark 
matter halos. We vary the redshift, width and size of the drop in the SFR. We 
assume that 50% of the high-redshift SNe are Type IIP and 50% are Type IIL. 
We make use of the observed properties of local SNe, such as their lightcurves 



(|26l ). spectra (|27t). and peak magnitude distributions (28), to predict the cor- 
responding number of detectable SNe as a function of redshift. Specifically, 
we consider JWST, a 6m diameter space telescope, scheduled for launch in 
2013. The relevant instrument on JWST is NIRcam, a near-infrared imaging 
detector with a FOV of 2.3' x 4.6' and simultaneous exposures in two filters. 

We find that 4-24 SNe may be detectable from z ^ 5 at the sensitivity of 
3 nJy (requiring 10 5 s exposures in a 4.5 fim band) in each ~ 10 arcmin 2 
JWST field. In a hypothetical one year survey (or a collection of ~ hundred 
repeated exposures), we expect to detect up to thousands of SNe per unit 
redshift at z ~ 6. Our results imply that, for most scenarios, high-redshift 
SNe observations can be used to detect reionization features out to z ~ 13, as 
well as set constraints on the photoionization heating feedback on low-mass 
halos at the reionization epoch. Specifically, for a wide range of scenarios at 
z TC ^ 13, the drop in the SNR due to a sharp feature in the reionization history 
can be detected at S/N ^ 3 with only tens of deep J WST exposures. In cases 
with strong negative feedback, the drop in the SNR can be detected if it's 
spread out over Az re ^ 4; less optimistic scenarios require Az re ~ 1 - 2. Even 
in the case of a more extended reionzation history, fairly sharp features are 
likely, which could leave detectable signatures in the SNR (see Fig. 9 in fill), 
and associated discussion). Our results therefore suggest that future searches 
for high-2 SNe could be a valuable new tool, complementing other techniques, 
to study the process of reionization, as well as the feedback mechanism that 
regulates it. 



5 Conclusions 



We have developed several useful tools designed to probe reionization. The 
main implications of these works can be summarized as follows. 

• The cosmological proximity effect can be used to study imprints of the 
host environment properties on the source's spectra; specifically, a neutral 
fraction of x^i = 1 and be distinguished from £ H i ^ 0.01 with only tens of 
bright QSO spectra (or ~ hundreds of spectra which are 100 times fainter), 
in the presence of emission uncertainties of ~ 20%. 

• The spectral analysis of QSO SDSS J1030+0524 suggests that the Universe 
(or at least a large part of it) is very neutral at z ~ 6.2, with x H i ^ 0.2. 
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• The ionizing luminosity of SDSS J1030+0524 is a factor of ~ 2-10 lower than 
expected, strengthening arguments against quasar-dominated reionization. 

• Disparate Lyman line absorption statistics can probe ionization topology. 

• The sharpness of the Stromgren sphere boundary can be used to constrain 
the X-ray contribution to the quasar's intrinsic emission. 

• Future instruments, such as JWST, can yield many high-z SNe; specifically, 
4-24 SNe/field may be detectable from z ^ 5 at the sensitivity of 3 nJy 

• High-z SNe can be used to detect reionization features out to z & 13. 

I thank my advisor, Zoltan Haiman, and my collaborators Renyue Cen and 
Benjamin Johnson for permission to draw on joint work. 



References 

[I] X. Fan, et al, AJ 123, 1247 (2002). 

[2] C. L. Bennett, et al, Astrophys. J. Suppl. 148, 1 (2003). 

[3] Z. Haiman and G. .P. Holder, Astrophys. J. 595, 1 (2003). 

[4] D. N. Spergel, et al, Astrophys. J. Suppl. 148, 175 (2003). 

[5] A. Mesinger, Z. Haiman, and R. Cen, Astrophys. J. 613, 23 (2004). 

[6] A. Mesinger and Z. Haiman, Astrophys. J. 611, 69 (2004). 

[7] R. L. White, et al, AJ 126, 1 (2003). 

[8] R. Barkana and A. Loeb, Astrophys. J. 601, 64 (2004). 

[9] P. Madau and M. J. Rees, Astrophys. J. 542, 69 (2000). 

[10] R. Cen and Z. Haiman, Astrophys. J. 542, 75 (2000). 

[II] M. Elvis, et. al, Astrophys. J. Suppl. 95, 1 (1994). 

[12] R. C. Telfer, W. Zheng, and G. A. Kriss, A J 565, 773 (2002). 

[13] M. Dijkstra, Z. Haiman, and A. Loeb, Astrophys. J. 613, 646 (2004). 

[14] J. S. B. Wyithe, and A. Loeb, Nature 427, 815 (2004). 

[15] J. A. Tyson, Proc. SPIE Int. Soc. Opt. Eng. 4836, 10 (2003). 

[16] A. Mesinger, B. D. Johnson, and Z. Haiman, preprint 

|arXiv:astro-ph/0505110| . 
[17] P. Madau et al, MNRAS 283, 1388 (1996). 
[18] R. Barkana and A. Loeb, Astrophys. J. 539, 20 (2000). 
[19] R. Cen and P. McDonald, Astrophys. J. 570, 457 (2002). 
[20] A. Mesinger, R. Perna, and Z. Haiman, Astrophys. J. 623, 1 (2005). 
[21] P. Madau, M. della Valle, and N. Panagia, MNRAS 297, 17 (1998). 
[22] J. H. Wise and T. Abel, preprint |arXiv:astro-ph/0411558 . 
[23] J. Miralda-Escude and M. J. Rees, Astrophys. J. 478, 57 (1997). 
[24] T. Dahlen and C. Fransson, AAP 350, 349 (1999). 
[25] A. A. Thoul and D. H. Weinberg, Astrophys. J. 465, 608 (1996). 
[26] J. B. Doggett and D. Branch, AJ 90, 2303 (1985). 
[27] E. Baron, et al, Astrophys. J. 545, 444 (2000). 
[28] D. Richardson, et al, A J 123, 745 (2002). 



10 



